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Gamma-ray bursts (GRBs) are among the most violent occurrences in the universe. They are 
powerful explosions, visible to high redshift, and thought to be the signature of black hole birth. 
They are highly luminous events and provide excellent probes of the distant universe. GRB research 
has greatly advanced over the past 10 years with the results from Swift, Fermi and an active follow- 
up community. In this review we survey the interplay between these recent observations and the 
theoretical models of the prompt GRB emission and the subsequent afterglows. 
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I. INTRODUCTION 

GRBs are the most extreme explosive events in the 
Universe. The initial prompt phase lasts typically less 
than 100 s and has an energy content of ~ 10 51 ergs, 
giving a luminosity that is a million times larger than 
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the peak electromagnetic luminosity of a supernova. The 
GRB name is a good one because their spectra peak in 
the gamma-ray band between ~ 100 keV and ~ 1 MeV. 
The source of the energy powering the bursts is thought 
to be the gravitational collapse of matter to form a black 
hole or other compact object. 

GRBs were discovered in the late 1960s by the Vela 
satellites monitoring the Nuclear Test Ban Treaty be- 
tween the US and the Soviet Union [1 . It was slow 
progress for 20 years learning the origin of these brilliant 
flashes. Gamma-ray instruments of the time had poor 
positioning capability so only wide-field, insensitive tele- 
scopes could follow-up the bursts to look for counterparts 
at other wavelengths. Nothing associated with the GRBs 
was seen in searches hours to days after their occurrence. 
In the 1990's the Burst and Transient Source Experiment 
(BATSE) onboard the Compton Gamma Ray Observa- 
tion (CGRO) obtained the positions of ~ 3000 GRBs and 
showed that they were uniformly distributed on the sky 
[2] , indicating either an extragalactic or galactic-halo ori- 
gin. BATSE also found that GRBs separate into two du- 
ration classes, short and long GRBs, with dividing line at 
~ 2 s [3 j . The detection of X-ray afterglows and more ac- 
curate localizations delivered by the BeppoSAX mission 
[4j [5] enabled optical redshifts of GRBs to be measured 
and their extragalactic origin to be confirmed. The long 
bursts were found to be associated with far-away galaxies 
at typical redshifts z = 1 — 2, implying energy releases 
in excess of 10 50 ergs, after taking into account the jet 
beaming factor. The afterglow time decay often steep- 
ened after a day indicating a geometrical beaming of the 
radiation into a jet of opening angle ~ 5° [6]. 

Key open issues prior to the Swift and Fermi era in- 
cluded: (1) the origin of short GRBs, (2) the nature of 
the high energy radiation from GRBs, (3) the redshift 
distribution of bursts and their usage for early universe 
studies, and (4) the physics of the jetted outflows. 

Swift and Fermi, launched in 2004 and 2008 respec- 
tively, have opened a new era in GRB research. 

Swift is a NASA mission and Fermi a NASA collab- 
oration with the US Department of Energy. They are 
both major international partnerships and have different 
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and complementary capabilities. Swift has a wide-field 
imaging camera in the hard X-ray band that detects the 
bursts at a rate of ~ 90 per year, providing positions with 
arcminute accuracy. The spacecraft then autonomously 
and rapidly (100 s) reorients itself for sensitive X-ray 
and UV/optical observations of the afterglow. Fermi has 
two wide-field instruments. One detects bursts in the 
gamma-ray band at a rate of ~ 300 per year, provid- 
ing spectroscopy and positions with 10-degree accuracy. 
The other observes bursts in the largely-unexplored high- 
energy gamma-ray band at a rate of ~ 10 per year. Com- 
bined, the two missions are advancing our understand- 
ing of all aspects of GRBs, including the origin of short 
bursts, the nature of bursts coming from the explosion of 
early stars in the universe and the physics of the fireball 
outflows that produce the gamma-ray emission. 



II. OBSERVATIONS 
A. Swift GRBs 

The Swift mission [7] has three instruments: the Burst 
Alert Telescope (BAT; 0), the X-Ray Telescope (XRT; 
[9]) and the UV Optical Telescope (UVOT; [TO]). The 
BAT detects bursts and locates them to ~ 2 arcminute 
accuracy. The position is then sent to the spacecraft 
to repoint the XRT and UVOT at the event. Positions 
are also rapidly sent to the ground so that ground tele- 
scopes can follow the afterglows. There are more than 
50 such telescopes of all sizes that participate in these 
world-wide follow-up campaigns. Measurements of the 
redshift and studies of host galaxies are typically done 
with large ground-based telescopes, which receive im- 
mediate alerts from the spacecraft when GRBs are de- 
tected. Swift has, by far, detected the largest number 
of well-localized bursts with afterglow observations and 
redshift determinations. As of 1 September 2012, BAT 
has detected 709 GRBs (annual average rate of ~ 90 per 
year). Approximately 80% of the BAT-detected GRBs 
have rapid repointings (the remaining 20% have space- 
craft constraints that prevent rapid slewing). Of those, 
virtually all long bursts observed promptly have detected 
X-ray afterglow. Short bursts are more likely to have neg- 
ligible X-ray afterglow, fading rapidly below the XRT 
sensitivity limit. The fraction of rapid-pointing GRBs 
that have UVOT detection is ~ 35%. Combined with 
ground-based optical observations, about ~ 60% of Swift 
GRB have optical afterglow detection. There are 225 
Swift GRBs with redshifts compared with 41 pre-Swift. 

Swift was specifically designed to investigate GRB af- 
terglows by filling the temporal gap between observations 
of the prompt emission and the afterglow [TT] . Figure [I] 
shows example X-ray afterglows from XRT from long and 
short GRBs [12]. The combined power of the BAT and 
XRT has revealed that in long GRBs the non-thermal 
prompt X-ray emission smoothly transitions into the de- 
caying afterglow. Often, a steep-to-shallow transition is 
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FIG. 1: Representative examples of X-ray afterglows of (a) 
long and (b) short Swift events with steep-to-shallow transi- 
tions (GRB050315, 050724), large X-ray flares (GRB050502B, 
050724), and rapidly declining (GRB051210) and gradually 
declining (GRB051221a, 050826; flux scale divided by 100 
for clarity) afterglows. From Gehrels, Ramirez- Ruiz and Fox 
(2009) Q2]. 



found suggesting that prompt emission and the after- 
glow are distinct emission components [13] [14]. The early 
steep-decay phase seen in the majority of GRBs is a sur- 
prise. The current best explanation is that we are seeing 
high-latitude emission due to termination of central en- 
gine activity [T4] . 

Swift has discovered erratic flaring behavior (Figure[T]), 
lasting long after the prompt phase in ~ 30% of X-ray 
afterglows. The most extreme examples are flares with 
integrated power similar to or exceeding the initial burst 
[T5] . The rapid rise and decay, multiple flares in the same 
burst, and cases of fluence comparable with the prompt 
emission, suggest that these flares are due to the central 
engine. 

Long GRBs (LGRBs) are associated with the bright- 
est regions of galaxies where the most massive stars oc- 
cur [IB] . LGRBs occur over a large redshift range from 
z = 0.0085 (GRB 980425) to z > 8 (GRB 090423 & GRB 
090429B). With few exceptions, LGRBs that occur near 
enough for supernova detection have accompanying Type 
lb or Ic supernovae. The exceptions are cases of bursts 
that may be misclassified as long or that may have ex- 
ceptionally weak supernovae. These facts support the 
growing evidence that long bursts are caused by "collap- 
sars" where the central core of a massive star collapses 
to a compact object such as a black hole [17] or possibly 
a magnetar [18] [19] . 

Theories for the origin of LGRBs predate Swift, but 
are supported by the new data. In them, a solar rest 
mass worth of gravitational energy is released in a very 
short time (seconds or less) in a small region of the order 
of tens of kilometers by a cataclysmic stellar event. The 
energy source is the collapse of the core of a massive 
star. Only a small fraction of this energy is converted 
into electromagnetic radiation, through the dissipation 
of the kinetic energy of a collimated relativistic outflow, 
a fireball with bulk Lorentz factors of T ~ 300, expanding 
out from the central engine powered by the gravitational 
accretion of surrounding matter into the collapsed core 
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or black hole. 



B. Fermi GRBs 

Mission & Statistics — The Fermi instruments are the 
Gamma-ray Burst Monitor (GBM, [20]) and the Large 
Area Telescope (LAT, [21]). The GBM has scintillation 
detectors and covers the energy range from 8 keV to 40 
MeV. It measures spectra of GRBs and determines their 
position to ~ 5° accuracy. The LAT is a pair conversion 
telescope covering the energy range from 20 MeV to 300 
GeV. It measures spectra of sources and positions them 
to an accuracy of < 1°. The GBM detects GRBs at 
a rate of ~ 300 per year, of which on average 20% are 
short bursts. The LAT detects bursts at a rate of ~ 10 
per year. 

LAT bursts have shown two common and interesting 
features: (1) delayed emission compared to lower energy 
bands and (2) longer lasting prompt emission compared 
to lower every bands. The four brightest LAT bursts are 
GRB 080916C [22], GRB 090510 [23[|24], GRB 090902B 
[25], and GRB 090926A [26 . They have yielded hun- 
dreds of > 100 MeV photons each, and together with 
the lower energy GBM observations, have given unprece- 
dented broad-band spectra. In GRB 080916C, the GeV 
emission appears only in a second pulse, delayed by ~ 4 
s relative to the first pulse (Figure |2|. Such a delay is 
present also in short bursts, such as GRB 090510, where 
it is a fraction of a second. This soft-to-hard spectral 
evolution is clearly seen in all four of these bright LAT 
bursts, and to various degrees a similar behavior is seen 
in other weaker LAT bursts. 

In some bursts, such as GRB 080916C and several oth- 
ers, the broad-band gamma-ray spectra consist of a sim- 
ple Band- type broken power-law function in all time bins. 
In GRB 0809 16C the first pulse has a soft high energy in- 
dex disappearing at GeV energies, while the second and 
subsequent pulses have harder high energy indices reach- 
ing into the multi-GeV range. In some other bursts, such 
as GRB 090510 [231 EU and GRB 090902B [25], a sec- 
ond hard spectral component extending above 10 GeV 
without any obvious break appears in addition to com- 
mon Band spectral component dominant in the lower 8 
keV-10 MeV band. 

An exciting discovery, unanticipated by results 
from the Energetic Gamma-Ray Experiment Telescope 
(EGRET) instrument on CGRO, was the detection 
of high-energy emission from two short bursts (GRB 
081024B [27] and GRB 090510 [23]|24]). Their general 
behavior (including a GeV delay) is qualitatively similar 
to that of long bursts. The ratio of detection rates of 
short to long GRBs by LAT is - 7% (2 out of 27), which 
is significantly smaller than the ~ 20% by GBM. 

While the statistics on short GRBs are too small to 
draw firm conclusions, so far, the ratio of the LAT flu- 
ence to the GBM fluence is ~ 100% for the short bursts as 
compared to ~ 5-60% for the long bursts. It is also note- 
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FIG. 2: Light curves of GRB 080916C with the GBM (top 
three panels) and LAT (bottom two panels) . The high energy 
LAT emission is delayed relative to the lower energy GBM 
emission. From Abdo et al. 2009 [22] . 



worthy that, for both long and short GRBs, the > 100 
MeV emission lasts longer than the GBM emission in the 
< 1 MeV range. The flux of the long-lived LAT emission 
decays as a power law with time, which is more reminis- 
cent of the smooth temporal decay of the afterglow X- 
ray and optical fluxes rather than the variable temporal 
structure in the prompt keV-MeV flux. This similarity 
in the smooth temporal evolution of the fluxes in differ- 
ent wave bands has been detected most clearly in GRB 
090510 [28, 29 although this burst also requires a sepa- 
rate prompt component to the LAT emission [30]. This 
short burst was at z = 0.9, and was jointly observed by 
the LAT, GBM, BAT, XRT and UVOT. 



The LAT detects only ~ 10% of the bursts triggered 
by the GBM which were in the common GBM-LAT field 
of view. This may be related to the fact that the LAT- 
detected GRBs, both long and short, are generally among 
the highest fluence bursts, as well as being among the 
intrinsically most energetic GRBs. For instance, GRB 
080916C was at z = 4.35 and had an isotropic-equivalent 
energy of E{ so « 8.8 x 10 54 ergs in gamma rays, the largest 
ever measured from any burst [22] . The long LAT bursts 
GRB 090902B [25 at z = 1.82 had E iso « 3.6 x 10 54 ergs, 
while GRB 090926A [26 at z = 2.10 had E iso « 2.24 x 
10 54 ergs. Even the short burst GRB 090510 at z = 0.903 



produced, within the first 2 s, an E^ 



1.1 x 10 53 ergs 
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C. High redshift GRBs 

Swift is detecting GRBs at higher redshift than pre- 
vious missions due to its higher sensitivity and rapid af- 
terglow observations. The average redshift for the Swift 
GRBs is 2.3 compared to 1.2 for previous observations. 
Although statistics are poor, the highest redshift GRBs 
are seen to have high luminosity, resulting in fluxes well 
above the detection threshold. Such bursts are also 
strong at other wavelengths. Table 1 presents optical 
data for the highest redshift GRBs observed to date. The 
optical brightness is given in magnitudes in R, /, J and 
K wavelength bands. The bursts are 3-4 orders of mag- 
nitude brighter than typical galaxies at these redshift. 

To date, Swift has seen seven GRBs with z > 5. 
The highest z GRBs so far are z = 6.3 (GRB 050904 
[31]), z = 6.7 (GRB 080913 [32]), z = 8.2 (GRB 090423 
[33] EU), and z = 9.4 (GRB 090429B [35]). Figure [| 
shows a spectrum of the current record holder with a 
spectroscopic redshift, GRB 090423 [33]. Because GRBs 
are so bright across the electromagnetic spectrum, they 
are excellent tools for studying the high redshift universe. 
The information that can be provided by GRBs is com- 
plementary to what can be learned by studies of galaxies 
and quasars. For example, at z > 8 it is expected that 
the bulk of the star formation activity is taking place in 
small galaxies below the detection limits of HST. 

The chemical evolution of the universe can be studied 
with GRBs as illustrated in Figure [4] [36 . GRBs provide 
data to higher redshift than active galaxies. The metal- 
licity bias of GRBs is currently contradictory. Plots like 
Figure [4] show GRBs in higher metallicity star forming 
regions of galaxies. However, comparing star forming 
regions, GRBs actually tend to favor ones with lower 
metallicity [37] . 

Another way that GRBs are contributing to our under- 
standing of the high-redshift universe is in the determi- 
nation of the star formation history (Figure [5]). LGRBs 
are the endpoints of the lives of massive stars and their 
rate is therefore approximately proportional to the star 
formation rate. They give information at high redshift 
where the rate is highly uncertain. There may be evo- 
lutionary biases, such as a dependence of LGRBs on the 
metallicity of host galaxies, so studies relating them to 
the star formation rate must include these factors into 
account [381 [39] . 



TABLE I: Optical Observations of High Redshift GRBs 
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FIG. 3: The composite spectrum of the GRB 090423 af- 
terglow obtained with the VLT. Also plotted are the sky- 
subtracted photometric data points obtained using Gemini 
North & South and VLT. A model spectrum showing the H 
I damping wing for a host galaxy with a hydrogen column 
density of NH I 5 x 10 21 cm -2 at a redshift of z = 8.23 is also 
plotted (solid black line). Inset, allowing for a wider range in 
possible host NH I values gives the la (68%) and 2a confi- 
dence contours shown. From Tanvir et al. 2009 33 . 
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FIG. 4: Redshift evolution of the metallicity (represented 
here by the ratio of oxygen to hydrogen abundance) relative to 
solar values, for GRBs shown with blue dots and active galaxy 
quasars show with open circles. The abundances are deter- 
mined from spectral absorption lines in the continuum radia- 
tion. GRB lines are predominantly from the gas in the host 
galaxy in the star-forming region near the explosion, whereas 
quasars are of random lines-of-sight through the galaxy. The 
GRB metallicity is on average ~ 5 times larger than in QSO. 
These are based on damped Lyman alpha (DLA) spectral 
features. The upper horizontal x-axis indicates the age of the 
Universe (Hubble time). From Savaglio et al. 2006 |36| . 
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FIG. 5: Cosmic star formation history. Shown are the data 
compiled in Hopkins & Beacom 40 (light circles) and con- 
tributions from Ly emitters (LAE) [41]. Recent LBG data 
are shown for two UV LF integrations: down to 0.2L* (down 
triangles; as given in Bouwens et al. [42]) and complete z = 3 
(up triangles). Swift GRB-inferred rates are diamonds, with 
the shaded band showing the range of values resulting from 
varying the evolutionary parameters. Also shown is the criti- 
cal density derivative from Madau et al. [43] for C / f esc = 40, 
30, 20 (dashed lines, top to bottom). From Kistler et al. 2009 
|38j . See also study by Robertson & Ellis [39] . 



D. GRB-Supernova connection 

On 18 February 2006 Swift detected the remarkable 
burst GRB 060218 that provided considerable new infor- 
mation on the connection between SNe and GRBs. It was 
longer (35 min) and softer than any previous burst, and 
was associated with SN 2006aj at only z = 0.033. SN 
2006aj was a (core-collapse) SN Ib/c with an isotropic 
energy equivalent of a few 10 49 erg, thus under lumi- 
nous compared to the overall energy distribution for long 
GRBs. The spectral peak in prompt emission at ~ 5 keV 
places GRB 060218 in the X-ray flash category of GRBs 
[39] , the first such association for a GRB-SN event. Com- 
bined BAT-XRT-UVOT observations provided the first 
direct observation of shock-breakout in a SN [44] . This is 
inferred from the evolution of a soft thermal component 
in the X- ray and UV spectra, and early-time luminosity 
variations. Concerning the SN, SN 2006aj was dimmer 
by a factor ~ 2 than the previous SNe associated with 
GRBs, but still ~ 2-3 times brighter than normal SN Ic 
not associated with GRBs [45] [46]. GRB 060218 was an 
underluminous burst, as were two of the other three pre- 
vious cases. Because of the low luminosity, these events 
are only detected when nearby and are therefore rare. 
However, they are actually ~ 5-10 times more common 
in the universe than normal GRBs [47] and form a dis- 
tinct component in the luminosity function [48] . 

GRB 060218 joins 4 other previous cases of GRBs from 
other missions that were associated with specific super- 



novae (980425, 021211, 030329, 031203). In addition, 
there are several other Swift cases (050525 A, 091127, 
100316D, 101219B and 120422A). Table 2 gives a sum- 
mary of all long GRBs that have had either specific super- 
novae or features in their afterglow light curves suggestive 
of the presence of a supernova. Many of the nearby GRBs 
with specific supernovae are in the underluminous class 
like GRB 060218. There are also a few cases listed in the 
table of nearby GRBs that had supernova searches, with- 
out a detection. Two of these, GRB 060505 and 060614 
are discussed further in section Hi Fl 

TABLE II: Nearby GRBs and Supernova Detections or Limits 



GRB 


Redshift, 


z Type a 


SN Search 


980425 


0.0085 


long-UL 


oat -i AAoi 
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020903 


0.251 


XRF-UL 


LC bump & spectrum 
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0.55 
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0.283 


long-UL 


NS 2012bz 
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0.36 


long-UL 


LC bump & spectrum 
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0.297 


long-UL 


LC bump 
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0.54 


long- UL 


LC bump 
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0.54 


long 


LC bump 
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LC bump GROND 
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long 


LC bump GROND 


090424 




long 


LC bump GROND 


100902A 




long 


LC bump GROND 


110402A 




long 


LC bump GROND 


040701 


0.215 


XRF-UL 


no SN (< 0.1 SN98bw) 


060505 


0.089 


"long"-UL 


no SN (< 0.004 SN98bw) 


060614 


0.125 


"long" 


no SN (< 0.01 SN 98bw) 


101225A 


0.40 


long 


no SN (GCN 11522) 


a UL = underluminous, XRF = X-Ray Flash 



E. Short GRBs 

At the time of Swifts launch, the greatest mystery of 
GRB astronomy was the nature of short-duration, hard- 
spectrum bursts (SGRBs). Although more than 50 long 
GRBs had afterglow detections, no afterglow had been 
found for any short burst. In summer 2005 Swift and 
HETE-2, precisely located three short bursts for which 
afterglow observations were obtained leading to a break- 
through in our understanding of short bursts [49H55] . As 
of 2012, BAT has detected 65 SGRBs, 80% of which have 
XRT detections, and 15 of which have redshifts. 

In contrast to long bursts, the evidence is that SGRBs 
typically originate in host galaxies with a wide range of 
star formation properties, including low formation rate. 
Their host properties are substantially different than 
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those of long bursts [56] [57] indicating a different origin. 
Also, nearby SGRBs show no evidence for simultaneous 
supernovae ([58]; and references therein), very different 
than long bursts. Taken together, these results support 
the interpretation that SGRBs arise from an old popula- 
tions of stars and are due to mergers of compact binaries 
(i.e., double neutron star or neutron star - black hole) 
[58H60] . It is possible that there is a subpopulation of 
high-redshift, high-luminosity short GRBs with a differ- 
ent origin [6T| 162]. 

Measurements or constraining limits on beaming from 
light curve break searches have been hard to come by 
with the typically weak afterglow of SGRBs. With large 
uncertainties associated with small number statistics, the 
distribution of beaming angles for SGRBs appears to 
range from ~ 5° to > 25° [65 , 66 , roughly consistent but 
perhaps somewhat larger than that of LGRBs. Swift ob- 
servations also revealed long (~ 100 s) "tails" with softer 
spectra than the first episode following the prompt emis- 
sion for about 25% of short bursts [67, 68 . Swift local- 
ization of a short GRB helps narrow the search window 
for gravitational waves from that GRB [69 . Detection of 
gravitational waves from a Swift GRB would lead to great 
scientific payoff for merger physics, progenitor types, and 
NS equations of state. 



F. Unusual GRBs and transients 

Here is a list of unusual transients detected by Swift 

Short GRB 050925: This unusual burst triggered the 
BAT with a single peaked outburst of duration T90 = 70 
ms [7TJ[72]. It occurred near the galactic plane, and noth- 
ing was seen in the UVOT. However, the F-band extinc- 
tion toward the source was AV = 7.05 mag. The XRT 
spectra and lightcurve show no significant X-ray emission 
in the field, suggesting that any X-ray counterpart to this 
burst was faint. Markwardt et al. (2005) [72] were able 
to fit a power law spectrum to the BAT data, with a pho- 
ton index 1.74 ±0.17; they found a better fit could be ob- 
tained with a blackbody spectrum kT = 15.4±1.5 keV, a 
value consistent with small- flare events from SGRs. The 
low galactic latitude and soft spectrum indicate a possi- 
ble galactic source or SGR. 

GRBs 060505 & 060614: GRB 060505 and GRB 
060614 were nearby GRBs (z = 0.089 and 0.125, resp.) 
with no coincident SN (Table 2), to deep limits [73]. GRB 
060614 was bright (15-150 keV fluence of 2.2 x l6~ 5 erg 
cm -2 ), and, with a T90 of 102 s, seemed to be a secure 
long GRB. Host galaxies were found [75H75] and deep 
searches were made for coincident SNe. All other well- 
observed nearby GRBs have had SNe, but GRB 060614 
did not to limits > 100 times fainter than previous de- 
tections [72H74] . GRB 060614 shares some characteris- 
tics with SGRBs [68]. The BAT light curve shows an 
initial short hard flare lasting ~ 5 s, followed by an ex- 
tended softer episode, ~ 100 s. The light curve is similar 
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FIG. 6: Spectral lag as a function of peak luminosity showing 
GRB 060614 in the region of short-duration GRBs [68]. 



to some Swift SGRBs and a subclass of BATSE SGRBs 
[67] . GRB 060614 also falls in the same region of the 
lag-luminosity diagram as SGRBs (Figure^. Thus GRB 
060614 is problematic to classify. It is a LGRB by the tra- 
ditional definition, but lacks an associated SN. It shares 
some similarities with SGRBs, but the soft episode is 
brighter, which would be difficult to account for in the 
NS-NS merger scenario. 

Hostless GRB 070125: There is no obvious host galaxy 
for GRB 070125. Deep ground-based imaging reveals no 
host to R > 25.4 mag. Cenko et al. (2008) [76] present 
an analysis of spectroscopic data which reveals weak Mg 
II lines indicative of halo gas. In the field are two blue 
galaxies offset by > 27 kpc at z = 1.55. If there is an 
association with one of them, it would imply a velocity 
~ 10 4 km s _1 over a ~ 20 Myr lifetime of the massive 
progenitor. The only known way of achieving this would 
have been a prior close interaction with a massive BH. 
However, this interpretation was muddied by Chandra et 
al. (2008) [77], who inferred a dense environment, based 
on bright, self-absorbed radio afterglow. They proposed 
a scenario in which the high density material lies close to 
the explosion site, and the lower density material further 
away. They note GRB 070125 was one of the brightest 
GRBs ever detected, with an isotropic release of 10 54 erg 
(by comparison, M c 2 ~ 2 x 10 54 erg). The prompt emis- 
sion from GRB 070125 was also seen by Suzaku/WAM 
[78]. 

Galactic GRB 070610: Discovered initially as GRB 
070610, this object, now dubbed Swift J195509.6+261406 
(Swl955+26), is thought to represent a member of a rel- 
atively new class of fast X-ray nova containing a BH. It 
had a duration of ~ 5 s, and also shows large variability. 
Kasliwal et al. (2008) [79 discuss several possibilities for 
this source, and propose an analogy with V4641 Sgr, an 
unusual BH binary which had a major outburst in 1999 
[80] . V4641 Sgr is a binary with a B9 III star orbiting a 
~ 9M@ BH [8J] and also exhibited strong and fast X-ray 
and optical variability. The analog is imperfect in that 
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the normal star in Swl955+26 is a cool dwarf rather than 
a B9 giant, suggesting a physical origin for the bursting 
behavior in the accretion disk and/or jet rather than the 
mass donor star. A ~ 5 s burst is certainly distinct from 
the ~ month long fast-rise exponential decay seen in sys- 
tems like A0620- 00 in 1975 [12183] which are thought to 
be due to a large scale storage and dumping of material 
in an accretion disk [84] [85], and may be more in line 
with either a disk-corona [86 or disk-jet instability [87] . 
Rea et al. (2011) [88] derive stringent upper limits on the 
quiescent X-ray emission from Swl955+26 using a ~ 63 
ks Chandra observation, and use this to argue against a 
magnetar interpretation. Simon et al. (2012) [89] present 
an analysis of optical emission during the 2007 outburst 
and find the optical emission manifests as spikes which 
decrease in duration as the burst progresses. They show 
that the emission can be explained by pure synchrotron 
emission. The main part of the optical outburst following 
the Swift/BAT trigger lasted ~ 0.3 d, with subsequent 
echo outbursts between 1 and 3 day, as in X-ray novae. 
The main part of the outburst has an exponential decay 
reminiscent of the 1975 outburst in A0620-00, but with a 
much shorter duration. If the global accretion disk limit 
cycle [84] [90] the mechanism for this outburst, it sug- 
gests a much shorter orbital period than the 7.75 hr for 
A0620-00 [91]. 

GRB 101225 Christmas burst: GRB 101225 was quite 
unusual: it had a T90 > 1700 s and exhibited a curv- 
ing decay when plotted in the traditional logF-logt co- 
ordinates. The total BAT fluence was ^ 3 x 10 -6 erg 
cm -2 . The XRT and UVOT found a bright, long-lasting 
counterpart. Ground based telescopes followed the event, 
mainly in R and /, and failed to detect any spectral fea- 
tures. At later times a color change from blue to red was 
seen; HST observations at 20 d found a very red object 
with no apparent host. Observations from the Spanish 
Gran Telescopio Canarias at 180 d detected GRB 101225 
at gAB = 27.21±0.27 mag and rAB = 26.90±0.14 [92]. 
Considered together, these characteristics are unique to 
this burst (Figure [7]), and led Campana et al. (2011) [93] 
to propose that it was caused by a minor body like an 
asteroid or comet becoming disrupted and accreted by a 
NS. Depending on its composition, the tidal disruption 
radius would be ~ 10 5 -10 6 km. Campana et al. find an 
adequate fit to the light curve by positing a ~ 5 x 10 20 
g asteroid with a periastron radius ~ 9000 km. If half 
the asteroid mass is accreted they derive a total fluence 
4 x 10 -5 erg cm -2 and distance ~ 3 kpc. Thone et al. 
(2011) [92] offer a different explanation for GRB 101225 
— the merger of a He star and a NS leading to a con- 
comitant SN. They derive a pseudo-redshift z = 0.33 by 
fitting the spectral-energy distribution and light curve 
of the optical emission with a GRB-SN template. Thus 
in their interpretation the event was much more distant 
and energetic. They argue for the presence of a faint, 
unresolved galaxy in deep optical observations, and fit 
the long-term light curve with a template of the broad- 
line type Ic SN 1998bw associated with GRB 980425. If 
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FIG. 7: Light curves of GRB 101225 in five energy bands: X- 
rays at 1 keV (black) , UV at 2030 angstroms (green) and 2634 
angstroms (blue), and optical at 6400 angstroms (R band, 
red) and 7700 angstroms (/ band, orange) [93]. The inset 
(with t in seconds) shows the Swift/XRT light curve. Shaded 
regions highlight the periastron passages calculated using a 
tidal disruption model [93 . 



their distance is correct as well as their interpretation 
of a component emerging at 10 d as being a SN, then 
its absolute peak magnitude MV a bs = —16.7 mag would 
make it the faintest SN associated with a long GRB. The 
isotropic-equivalent energy release at z = 0.33 would be 
> 1.4 x 10 51 erg which is typical of other long GRBs but 
greater than most other low-redshift GRBs associated 
with SNe. 



III. THEORY AND RADIATION MODELING 

Gamma-ray bursts are thought to be the endpoints of 
evolution of massive, > 3OM , stars (favorite for long 
GRBs) [12] or neutron star-neutron star/black hole close 
binaries (favorite for short GRBs) [58-60] [94]. In ei- 
ther case a black hole and a short-lived accretion disc 
are thought to be formed after the core collapse (long 
GRBs) or binary merger (short GRBs). In some cases a 
highly magnetized neutron star, called a magnetar, may 
form as an intermediate state which eventually collapses 
to a black hole [95]. Gravitational binding energy of the 
system is released within a short time scale which eventu- 
ally produces the GRB phenomenon when some fraction 
of the gravitational energy converts to gamma rays. This 
is done via accretion of material/gas onto the black hole. 

The compact core of a few solar mass with the 
Schwarzschild radius r g ~ 1O 6 (Mb#/3M ) cm forms the 
central engine of a GRB. Infalling gas onto the newly- 
formed black hole is accreted with high efficiency. Most 
of the gravitational energy, ~ 10 54 ergs, is radiated in 
~ 10 MeV thermal neutrinos as in core-collapse SNe and 
a lesser amount in gravitational wave. Only a small frac- 
tion, ~ 10 51 erg, is converted to a GRB fireball, which 
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eventually produces a highly-relativistic bi-polar jet and 
is responsible for the observed 7 rays [63] [64] . 

As discussed previously, observations of supernovae as- 
sociated with long GRBs provide strong support for the 
core-collapse schenario for this class of GRBs. The binary 
merger scenario of short GRBs has not been confirmed, 
but remains the most popular model. 

A. GRB fireball model 

The GRB fireball is produced when large amount of en- 
ergy is released within a short time scale. The isotropic- 
equivalent energy is E iso = (47r/Clj)Ej, where Ctj is the 
jet solid angle. The average solid angle (CI j) /Air ~ 1/500 
or the jet half-opening angle (0j) ~ 1/10 [Hll96] l97]. The 
luminosity of the outflowing material is L iso ~ E iso /t, 
where t is the typical duration of a long GRB. The tem- 
perature of the outflowing material at the base of the jet, 
with an inner radius of an accretion disc tq ~ 10 T cm, is 
kT = (L iso /ATirlca) ~ few MeV. Here a = 7.6 x lO" 15 
ergs cm -3 K -4 is the radiation density constant. The 
fireball contains gamma rays, electron-positron pairs and 
relatively fewer baryons (mostly protons). Due to the ra- 
diation pressure, the optically thick hot plasma expands 
and overcome the gravitational pressure and forms the 
relativistic jet [99] . 

An initially thermal energy dominated GRB jet be- 
comes kinetic energy dominated when the baryons carry 
the bulk of the initial fireball energy at a saturation ra- 
dius > ro [98 . The bulk Lorentz factor of the jet initially 
increases as T(r) ~ r/ro and becomes constant at a value 
r ~ r] no 10 2 -10 3 , where 77 ~ L iso /Mc 2 is defined as the 
ratio beteen the total energy to mass flow and is called 
the baryon loading of the fireball. The radius at which 
the fireball attains this coasting Lorentz factor is ~ tqT] 
[TOOl [TOT] . 

B. Prompt keV— MeV emission 

The prompt 7-rays from a GRB are thought to be 
emitted when a significant fraction of the fireball kinetic 
energy is converted back to radiation energy. The inter- 
nal shocks model is currently the most popular model 
for prompt 7-ray emission [102 . Collisions between 
shells of plasma or fireball produce shocks and a frac- 
tion of the bulk kinetic energy is converted into ran- 
dom particle energy which is then radiated, via syn- 
chrotron mechanism e.g., as the observed non-thermal 
7 rays. Given a timescale t v from observations of highly 
variable prompt 7-ray emission, internal shock collisions 
take place, due to a difference in Lorentz factor Ar ~ T 
between the shells, at a radius T{ ~ T 2 ct v ~ 3 x 
10 13 (r/316) 2 (^/0.01 s) cm. This is typically much larger 
than the photospheric radius of a baryon- dominated 
GRB jet, given by r ph ~ a T L iso /(4:7rT 3 m p c 3 ) ~ 4 x 
10 12 (L iso /10 52 erg s- 1 )(r/316) 3 cm (see, however, [T03] 



for different scaling of the photosphere). The observed 
pulse structure can be produced in colliding shells model 
|104[ 105 . Note that the variability time scale can be as 
small as ~ 1 ms, corresponding to the size scale of the 
jet base ro- 

The internal shocks are semi-relativistic, with Lorentz 
factor T s h ~ few. A fraction of the shock energy is 
thought to amplify magnetic field, via some kind of 
instabilities, up to an energy that is at some equipar- 
tition level cb with the jet luminosity Li so |106( 1107] . 
Another fraction of the shock energy is thought to ac- 
celerate electrons, via a Fermi mechanism, to relativistic 
energies. In a fast-cooling scenario, where electrons lose 
most of their energy to the observed radiation, e.g. via 
synchrotron/ Compton mechanism, within the dynamic 
time scale, the prompt 7-ray luminosity L 1 ^ so can be 
used to infer L iso (and eventually E iso from the duration 
of the burst) as L iso = L 1 ^ so /e e . The magnetic field 
in the jet frame can also be inferred from the prompt 
7-ray observations as B' = (2cbL 1 ^ so / e e r 2 Y 2 c) 1 ^ 2 ~ 2 x 

10 4 (e B /ee) 1/2 (L 7) Wl0 51 erg s- 1 ) 1 / 2 ^ /31Q)- 3 (t v /0-01 s)" 1 
Gauss |108] , With a minimum Lorentz factor of the 
electrons, 7 e ~ m p /m e , characteristic synchrotron 
radiation peak can explain the observed e p k ~ few 
hundred keV typical vF v peak of the 7-ray spectrum. 
The power-law part of the observed Band spectra above 
e p k can be explained as due to synchrotron radiation 
from the power-law distributed electron Lorentz factor, 
n( le )^ 1 -P^sn(e)^e-P/ 2 - 1 . 

The synchrotron shock model as described above has 
difficulties to explain the observed 7-ray spectrum below 
e p k in a large fraction of the bursts [10914111] . The fast- 
cooling synchrotron model is expected to produce a spec- 
trum n(e) no £~ 3 / 2 below e p k as contrary to e -1 , typical 
observed value. A harder, e -2 / 3 , synchrotron spectrum 
from electrons in random megnetic field is often called 
the "synchrotron death line" [109] which too has been 
observed violated in a sizeable fraction of GRBs. 

There have been many attempts to solve this Band 
spectrum crisis within the synchro-Compton mechanism. 
Explanations include synchrotron self-absorption in the 
keV range [112] , Compton scattering of optical photons 
to keV energies [113] . low-pitch angle scattering or jitter 
radiation [114(1115] , time-dependent acceleration [116] as 
well as mixing of keV photons with axion-like particles 
[117] , Photospheric thermal emission, as we decribe next, 
holds promise to explain the hard Band spectrum below 

C. Photospheric emission 

In case of a baryon-loaded GRB fireball, as we dis- 
cussed before, most of the thermal energy converts to 
kinetic energy at a radius r > rorj. The photospheric 
radius r p h is also typically much larger than 7*0 77. In case 
the fireball is highly baryon dominated r p h is smaller and 
could become comparable to tqT]. A significant fraction 
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of the thermal energy can be released from the fireball 
in this case [100 . Moreover dissipation of kinetic en- 
ergy at or below the photosphere, by shocks or other 
mechanisms such as neutron-proton decoupling and in- 
teractions |101[ I118H120| , could also enhance thermal ra- 
diation for moderately baryon-loaded fireballs as well as 
explain highly-efficient 7-ray emission from GRBs, which 
is difficult to explain with internal shocks. 

Simple photospheric models should produce thermal 
emission peaking at ~ 1 MeV, the temperature of the 
fireball at the base of the jet, which is in the range of ob- 
served Band peaks |121( 1122] . This leads to a plausible 
physical explanations for the Amati [123 or Ghirlanda 
|124j relations between the vF v peak energy and the 
total 7-ray energy |121[ 1125] . The temperature of the 
photosphere can be much below ~ 1 MeV in more com- 
plicated models, see e.g. [103 . Neutron-proton decou- 
pling, their interactions and resulting radiation from cas- 
cade particles at the photosphere or Compton scatter- 
ing of thermal photons by relativistic electrons in sub- 
photospheric shocks have also been explored to produce 
the phenomenological Band spectrum. These processes 
explain the power-law spectrum above s p k. The "syn- 
chrotron death line" issue [109 below e p k can be avoided 
in photospheric models. 

Recent modeling of Fermi data from GRB 090902B 
provides some evidence of photospheric emission being 
dominant in this burst [126] . In some other Fermi bursts 
a photospheric component in addition to the Band or 
Comptonized spectrum seems to improve fits, although 
thermal emission is not dominant in these bursts |127( 

[EH]. 



D. Poynting flux dominated GRB model 

Fermi observations of GRBs over a broad energy range 
and non-detection of a strong photospheric component 
in almost all bursts, and specially after GRB 080916C 
[2211129] , have renewed interests in magnetic or Poynting- 
flux dominated GRB jets, . A large dissipation radius, 
r ~ 10 15 cm, yet high efficiency inferred from modeling 
Fermi data can also be explained more naturally with 
magnetic models. Particles are accelerated in the recon- 
nection region of the magnetic field lines and produce 
the observed radiation. Magnetic models can be gen- 
eralized into two categories. In the first category the 
baryons are completely absent in the jet or dynamically 
negligible, at least initially |130H132] . and in the second 
category the baryon load is significant but dynamically 
sub-dominant relative to the magnetic stresses [133 -136 . 
A hybrid model, dubbed ICMART [T37] involves an ini- 
tially Poynting-flux dominated outflow with few baryons 
that lead to internal shocks when the magnetic energy 
subsides at a large radius. 



^ 1Q- 




FIG. 8: Prompt emission spectra of GRB 090926A at dif- 
ferent time intervals (bottom panel) as well as during the 
entire prompt phase (top panel). A power-law component in 
addition to the Band spectrum, as also found in other bright 
LAT-detected GRBs, is required to fit spectra. A cutoff in the 
PL component at few GeV has been detected in this GRB and 
in GRB 110731A [138]. From Abdo et al. [26]. 



E. GeV emission and plausible interpretations 

The origin of > 100 MeV radiation from GRBs and 
associated additional spectral component (Figure [8| has 
been intensely debated in recent years. In the context of 
the internal shock scenario, Synchrotron-Self-Compton 
(SSC) scattering of soft target photons by relativistic 
electrons is the first choice to produce the hard power- 
law component in a simple one-zone model [24 \ 1391 4142] . 
The observed time delay (Figure [2| between the keV- 
MeV radiation and > 100 MeV radiation, however, is 
longer than the expected delay for the SSC emission, 
which can be of the order of the variability time scale 
(< 1 s for long bursts) at best. Compton radiation from 
late internal shocks, with much longer variability time at 
GeV than in keV-MeV, can explain delayed LAT onset 
in some cases |140[ 1142] . 

An alternate mechanism to produce > 100 MeV ra- 
diation in the prompt GRB phase is through hadronic 
emission, either proton-synchrotron radiation and radia- 
tion from 77 — > e + e~ pair cascades [143] (Figure [9| or 
photohadronic interactions and associated cascade radi- 
ation |144| 1145] (Figure 10). In the one-zone radiating 
shells, the hadronic models can account for the delayed 
detection of the hard power-law component through re- 
quired proton acceleration and cooling time. 

The total energy of the burst, however, needs to be 
rather large in hadronic models because of a high jet 
Lorentz factor, T > 500, required to avoid 77 pair 
production in-situ by > 1 GeV photons. In particu- 
lar, to generate the observed 2.44 s delay in LAT for 
GRB 11 0731 A the comoving frame jet magnetic field 
needs to be B' « 10 5 (r/500)- 1/3 (t onset /2.44 s)~ 2 / 3 
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FIG. 9: Proton-synchrotron radiation model for delayed on- 
set of LAT emission in GRB 080916C. The delay is caused by a 
build-up of proton-synchrotron flux and a shift of cooling fre- 
quency from very high eenrgies to the LAT range, with time. 
GBM radiationn is assumed to be dominated by electron- 
synchrotron radiation. From 
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FIG. 10: Photohadronic cascade radiation model for the ad- 
ditional power-law component detected in GRB 090510. Syn- 
chrotron radiation from charged pions, muons and e + e~ con- 
tribute to the observed spectra. High-energy neutrinos are 
also expected in this model. From [144] . 



(Ky/100 MeV) x / 3 G, in case of proton-synchrotron 
model. The corresponding isotropic-equivalent jet lu- 
minosity is L jet > 10 57 (r/500) 16 / 3 (t onset /2.44 s) 2 / 3 
(£ 7 /100 MeV)" 2 / 3 erg s" 1 [131 CLE]. While the 
hadronic models require a larger total energy than the 
leptonic models, the energy budget can be brought down 
for T < 500, see also [T46l[T47] . 

A high T is not a prerequisite to avoid 77 absorption in 
case of two-zone radiation, in either leptonic or hadronic 
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FIG. 11: Light curves of > 100 MeV radiation from a few 
LAT-detected GRBs, which show decaying LAT flux as a com- 
mon feature. The grey band with slope £ -10 / 7 is the expected 
decay behavior for afterglow emission from a radiative blast 
wave. From [151]. 



models, where LAT emission comes from late internal 
shocks [HI [1431 ttH [US- The late internal shocks, how- 
ever, can not explain temporally extended LAT emission 
which is apparently non-variabile. 

The delayed onset, longer duration and smooth tem- 
poral decay of the > 100 MeV flux (Figure 11) led to 



the development of an early afterglow interpretation of 
high-energy radiation from the external forward shock 
[H [29j H50HT52] . The time delay between keV-MeV 
radiation and > 100 MeV radiation is explained as the 
time required for deceleration of the GRB jet with high T. 
The hard power-law component and temporally-extended 
LAT emission are interpretated as synchrotron radiation 
from the external shock while keV-MeV emission from 
the internal shocks in these two-zone emission models. 
Detailed modeling, however, shows that an extrapolation 
of the external shock emission at earlier time fails to re- 
produce all of the > 100 MeV flux in the prompt phase, 
requiring a significant contribution from the internal dis- 
sipation regions [30l I153H155] . 

Both in the internal- and external-shock scenario, mod- 
eling of high-energy radiation by the leptonic and/or 
hadronic processes depends crucially on T. Currently 
77 — » e + e~ pair-production constraint in internal shocks 
is the most widely used estimate of a lower limit on 
T [SHI [1561 IT57] . In the case of GRB 090926A and 
GRB 11 0731 A a cut off in the additional power-law com- 
ponent allowed to estimate T rather than its lower limit, 
assuming that the cut off at GeV energy is due to at- 
tenuation of flux by 77 — ^ e + e~ pair absorption. Figure 
[l2l shows the estimates of T from Fermi LAT-detected 
bright GRBs. There are, however, large (a factor 2-3) 
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FIG. 12: Bulk Lrentz factors of Fermi-LAT detected GRB 
jets calculated from 77 — > e + e~ pair absorption effect of GeV 
photons assuming that they are produced in the same physical 
region as the lower, < 100 MeV, photons. The lower limits, 
r m i n , are calculated in cases where there were no cutoff in 
the spectra in the GeV range. In the "PL only" lower limits 
On 1 rnin were calculated using only the additional polwer-law 
components of the spectra. Estimates of Y rather than r m i n 
were made when the observed > 1 GeV cutoff in the spectra 
were asumed due to 77 — >> e + e~ pair absorption. 
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FIG. 13: Multiwavelength light curves of GRB 110731A 
from Fermi GBM and LAT; Swift BAT, XRT and UVOT; 
and ground-based optical telescopes. Smooth decay of the 
fluxes across the wavelengths after ~ 6 s suggests their origin, 
including LAT emission, as from a decelerating blast wave in 
stellar wind environment. From 138 . 



dieted detection allowed the first measurements of host 
galaxies and redshift distances [4j [5] . 



uncertainty due to theoretical modeling of target radia- 
tion fields and environment [24j 115814160] . 



F. GRB afterglow 

Regardless of the details of the dissipation mecha- 
nism of the GRB jet to produce prompt 7-ray emis- 
sion, the accumulated surrounding material by the 
jet drives a blast wave and external forward shock 
[161] . A reverse shock is also expected in the baryon- 
loaded fireballs. The deceleration takes place with 
the jet energy compoarable to the kinetic energy of 
the swept-up material in the blast wave, Ei SO = 
(4/3)7rr^ ec n ext m p c 2 r 2 . The deceleration radius is r^ ec ~ 
2xl0 17 (iWl0 53 erg)V2( next /i cm^-V^r/lOO) 2 / 3 cm 
for a constant density external medium. 

After deceleration, the blast wave evolves in a self- 
similar fashion with the bulk Lorentz factor decreasing 
as T ^ t -3 / 8 for an adiabatic fieball and T ~ t~ 3 / 7 
for a radiative fireball, both in constant density medium 
[161, 162 . Particles are expected to be accelerated in 
the forward shock, which is highly relativistic. Magnetic 
field also assumed to be generated in the forward shock, 
which leads to synchrotron radiation by relativistic par- 
ticles |162[I163] . As the jet continues to be decelerated in 
the course of sweeping up more and more external mat- 
ter, the peak of the synchrotron flux shifts from high- 
frequency to low- frequency bands. Smooth decay of the 
flux across different energy bands is a characteristic fea- 
ture of the afterglow synchrotron model [164] , whose pre- 



Broadband X-ray to radio emission from GRBs has 
been quite successfully described by synchrotron radi- 
ation from a decelerating blast wave with or without 
variation of the straightforward models. Prompt opti- 
cal emission [164 , at the time the deceleration has been 
detected, with robotic ground-based telescopes such as 
ROTSE and others. The external shock interpretation 
of the late afterglow has proven robust overall, although 
some detailed features such as the X-ray steep decays fol- 
lowed by flat plateaus and occasional large flares are not 
well-understood yet. 



As discussed in Section |III E[ temporally extended 
> 100 MeV emission detected by Fermi-LAT from many 
GRBs well after the prompt keV-MeV emission has been 
advocated as afterglow emission by a number of authors 
[281 [291 H50HT52] . To test an early afterglow interpre- 
tation of the LAT data from GRBs one needs to con- 
front models with simultaneous multiwavelength data 
sets. Most recently, modeling of multiwavelength light 
curves (Figure fl3| and broadband Spectral Energy Dis- 
tribution (SEDjfrom GRB 110731A (Figure[l4| strongly 
indicate that [138 (i) LAT extended emission as early as 
To + 8.3 s, from broadband SED, is compatible with af- 
terglow synchrotron emission in other bands, (ii) A wind 
afterglow model [165-167 is favored over an ISM model, 
from the behavior of the light curves and SEDs. Detailed 
modeling of GRB 110731 A data also resulted in T ~ 500 
just before deceleration of the blast wave [138 . An in- 
ternal shock contribution to LAT data before To + 8.3 s 
can not be ruled out, however, and requires detailed in- 
vestigation with future Swift-Fermi bursts. 
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FIG. 14: Broadband spectra of GRB 110731A at two time 
intervals after the trigger. Single power laws fitting the spec- 
tra, together with flux decay behavior in Fig. [13] imply that 
LAT emission is afterglow. From [138] . 



G. High-energy neutrinos from GRBs 

Neutrinos with energies > 100 TeV are expected from 
photohadronic (jry) interactions of relativistic protons ac- 
celerated in the fireball internal or external shocks. In 
the case of internal shocks, the interaction is expected 
to take place between MeV photons produced by radi- 
ation from accelerated electrons and protons which are 
co-accelerated in the same shocks [168H171] . Neutrinos 
are expected to be produced more efficiently at an en- 
ergy e v > 100(£: 7 , p/c /100 keV)- 1 (r/316) 2 TeV. Neutrino 
production in internal shocks is expected to suppress at 
an energy e v > 1 PeV due to synchrotron and adiabatic 
losses of the high-energy pions and muon |169[ I172[ 1173] . 
In the case of external shocks, optical/UV photons pre- 
dominantly take part in pj interactions to produce > 1 
EeV neutrinos [TMnSJ. 

The most intense neutrino emission at < 10 TeV en- 
ergies from GRBs or collapsars might be due to 7 and 
pp interactions occurring inside the stellar envelope while 
the jet is still burrowing its way out of the star ]177| 1178] 
and before it successfully break through the stellar enve- 



lope to produce a GRB or choke inside. Figure 15 shows 
different neutrino flux models [179] from a nearby GRB 
030329 at z = 0.17 which was also connected with a su- 
pernova. 

The IceCube Collaboration has recently constrained 
high-energy neutrino emission from the internal shocks 
by stacking electromagnetic data from GRBs [180-182 . 
This constraint, in the ~ 1-20 PeV energy range, cal- 
culated from 40 and 59 string data severely excludes ef- 
ficient neutrino production in the internal shocks. This 
could be due to the high bulk Lorentz factor, > 500, 
of the GRB jet which reduces photopion production ef- 
ficiency, e.g. [147] 1183] . or due to inefficient aceleration 
of protons in the internal shocks. To accommodate Ice- 
Cube non-detection, revisions to the neutrino flux from 




6 7 8 

logl0[£ v /GeV] 

FIG. 15: Predicted neutrino fluxes from a relatively nearby 
GRB 030329 at a redshift z = 0.17. An accompanying super- 
nova was detected from this long GRB as well. From [179] , 



internal shocks have also been proposed [184] 1185] . New 
calculations in the models alternate to internal shocks 
are being done as well, e.g. [186 . Data from the full 
IceCube detector will either constrain the neutrino pro- 
duction efficiency in the internal shocks or rule out these 
models. Note that TeV neutrinos from the buried GRB 
jet [177(1178] and EeV neutrinos from the afterglow [174J - 
1176] remain viable. 



H. Gravitational waves from GRB 

Binary mergers of NS-NS and/or NS-BH, the most 
popular model for short GRBs, are the strongest can- 
didate sources of gravitational waves (GWs) with a total 
energy reaching ~ 10 54 ergs |187] 188 . The luminosity 
of GW from long GRBs, originating from core-collapse of 
massive stars, is highly model-dependent [189 . Recent 
numerical calculations of GW emission from long GRBs 
range from pessimistic [190] to modest [191] predictions. 
The expected detection rates of GW from compact bi- 
nary mergers in coincidences with short GRBs have been 
estimated to be several per year [192 for the advanced 
LIGO and VIRGO. Coordinated detection of GW with 
electromagnetic signal helps to reduce background signif- 
icantly, if short GRBs are the sources. Combining sub- 
threshold signal from electromagnetic, neutrino and GW 
observatories can also lead to significant detection 



IV. FUTURE GRB OBSERVATIONS 

Both Swift and Fermi are healthy and adequately- 
funded by NASA and foreign partners. Neither has life- 
time limits from expendable propulsion gasses or cryo- 
gens, and both have orbital lifetimes beyond 2025. They 
are joined by INTEGRAL (e.g., Mereghetti et al. 2003 
[194]) and MAXI [195], which also detect GRBs and the 
Interplanetary Network of GRB detectors [196 that lo- 
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calize events by triangulating their arrival times. 

There are several future missions with GRB capabili- 
ties that are being developed or proposed. They include: 

1. SVOM - A Chinese-French mission to detect GRBs 
and perform rapid optical and X-ray follow-up 

una 

2. UFFO - A Korean pay load flying on the Russian 
Lomonosov mission [198 . The payload includes a 
wide-field hard X-ray burst detector and a rapidly- 
orienting optical telescope. The pathfinder version 
of the mission will fly in 2012. 

3. Lobster - A wide-field focusing X-ray instrument 
operating in the 0.4-4 keV band. There are sev- 
eral versions that have been proposed including a 
UK instrument on the International Space Station 
[199 , a NASA Explorer mission including an after- 
glow IR telescope |200| , and an Explorer Mission 
of Opportunity on the International Space Station 
[201]. 

4. LOFT - A large-area X-ray timing mission selected 
for a Phase A study for ESA's M2 opportunity 
[202] , 

5. ASTROSAT - An Indian mission with hard X-ray 
imager, scanning sky monitoring, X-ray and UV 
telescopes [203 . 

Fermi LAT detection of > 10 GeV photons indicate 
that the prospects for detection of VHE radiation from 
GRBs by ground-based telescopes are promising. Highly 
sensitive Imaging Atmospheric Cherenkov Teleescopes 
(IACTs) such as MAGIC, VERITAS and HESS with a 
threshold energy as low as 25 GeV (special GRB mode 
for MAGIC) are currently slewing to the GRB direc- 
tions. Water Cherenkov detector HAWC is also un- 
der construction. Simultaneous observations by Fermi 
LAT and Cherenkov telescopes in the overlapping energy 
range will be extremely helpful to understand the nature 
of high-energy radiation from GRBs. Such observations 
can also pose strong constraints on both the leptonic and 
hadronic models, especially independent handles on the 
magnetic field strength are possible. Long-sought SSC 
radiation from the external forward shock or hadronic 
radiation in the prompt phase can also be searched for 
and/or constrained. Detection of very-high energy radi- 
ation from GRBs can be modeled to estimate the total 
explosion energy, which has implications for Graviational 
Wave detection by LIGO/VIRGO. 

There is a great opportunity for new science com- 
ing from future observations coordinated with GRBs of 
gravitational waves, neutrinos and across the electromag- 
netic spectrum. There are major new facilities coming 



on line and planned such as LIGO/VIRGO, IceCube, 
PanSTARRS, LSST and the 30m telescopes. 

V. CONCLUSIONS 

We are in an era of rapid progress in understanding 
GRBs. Swift, Fermi, INTEGRAL, MAXI, and IPN are 
all detecting GRBs. A vibrant community of observers is 
following up the GRB detections with new and capable 
instruments across the electromagnetic spectrum. Key 
findings from the past several years has included the fol- 
lowing: 

1. X-ray afterglow features, steep and shallow decays, 
flares, 

2. afterglows and accurate localizations of short bursts 
supporting a possible merger origin, 

3. confirmation of long GRB-Supernova connection, 

4. detection of high-redshift, z > 5, GRBs, 

5. delayed onset of GeV emission, 

6. high bulk Lorentz factor of the jet. 

With all this progress, there are many objectives that 
we hope to accomplish in the next 10 years. These in- 
clude: 

1. confirmation of short burst origin, 

2. jet opening angle and measurements of lighcurve 
breaks 

3. metallicity dependence 

4. relation to star formation history 

5. nature and interpretation of high energy spectral 
components. 

We look forward to many years more of GRB discovery. 
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